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Abstract: The strength of optomechanical interactions in a cavity optome- 
chanical system can be quantified by a vacuum coupling rate go analogous 
to cavity quantum electrodynamics. This single figure of merit removes 
the ambiguity in the frequently quoted coupling parameter defining the 
frequency shift for a given mechanical displacement, and the effective mass 
of the mechanical mode. Here we demonstrate and verify a straightforward 
experimental technique to derive the vacuum optomechanical coupling rate. 
It only requires applying a known frequency modulation of the employed 
electromagnetic probe field and knowledge of the mechanical oscillator's 
occupation. The method is experimentally verified for a micromechanical 
mode in a toroidal whispering-gallery-resonator and a nanomechanical 
oscillator coupled to a toroidal cavity via its near field. 
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1. Introduction 



Optomechanical systems are being explored for a variety of studies, ranging from exploring 
quantum limits of displacement measurements (Tpl, ground-state cooling of mechanical oscil- 
lators [3 ,4] to low-noise radiation-pressure driven oscillators J5}. In this new research field of 
cavity optomechanics |6j[7], rapid progress towards these goals has recently been enabled by 
the conception and realization of a large variety of nano- and micro-optomechanical systems, 
which implement parametric coupling of a mechanical oscillator to an optical or microwave 
resonator (for an overview of recently studied systems, see e. g. reference pi). Thereby, the 
resonance frequency co c of the electromagnetic cavity is changed upon a mechanical displace- 
ment x according to (0 c (x) = ft) c (0) +Gx, with the coupling parameter 

ax 

The resulting optomechanical interaction can be described by the interaction Hamiltonian 

0' t 

#int = fiGxh = PiGx^^ + a m ) ald D , (2) 

where x = x zp f(al n + a m ) is the mechanical displacement operator, h Q — aja the photon occu- 
pation of the optical cavity, and a , aj, a m and are the usual ladder operators of the optical 
(o) and mechanical (m) degrees of freedom. The root-mean-square amplitude of the mechanical 
oscillator's zero-point fluctuations 



Xzpf y 2m e ff£2 m ^ 

is determined by its resonance frequency Q. m and its effective mass m e ff as discussed in more 
detail below. 

The canonical example of such a system would be a Fabry-Perot resonator, one mirror of 
which is mounted as a mass on a spring moving as a whole along the symmetry axis of the cav- 
ity. In this case, the coupling parameter is given by G = —co c /L, where L is the cavity length, and 
the oscillator's effective mass simply is the moving mirror's mass. In most real optomechanical 
systems, however, internal mechanical modes constitute the mechanical degree of freedom of 
interest (for simplicity, our discussion will be restricted to a single mechanical mode only). In 
general, the displacement of such a mode has to be described by a three-dimensional vector 
field u(r). For every volume element located at r, it yields the displacement u(r) from its equi- 
librium position. In order to retain the simple one-dimensional description of equation Q, it is 
then necessary to map the vector field u{f) to a scalar displacement x. 

A frequently adapted procedure to address this mapping is to deliberately choose a coupling 
parameter G adequate to the overall geometry of the system, such as G = —(Q c /L for a Fabry- 



Perot-type cavity of length L (10 ■ 12 1, or G = —C0 C /R for a whispering-gallery-mode (WGM) 
resonator of major radius R ||8||13[. The equivalent displacement x is then defined by the meas- 
ured resonance frequency shift A© c induced by the displacement field u(r), i. e. x= Aco c /G. 
Following a different approach, the equivalent displacement x is defined, for example, as the 
displacement of the structure at a given position, x = |S(ro)| or the maximum displacement of 
a mode fl4Tfl7). Then , the coupling parameter G is calculated via the relation x = A(Q C /G. In 
some cases 1 14| 17 18, the coupling parameter is expressed as the effective length — i ^ 



of a resonator for which \G\ = <a c /L e ff- In this approach, the coupling parameter may be defined 
individually for each mechanical mode. 

In all cases, the only physically significant relation is the induced resonance frequency shift 
Ao) c upon a displacement u{f), which in a perturbation theory treatment, can be calculated 



according tcQ 

AcOc 1 J\E(r)\ 2 ■(e(r + u(r))-e(r))dr 3 



(4) 

Oh 2 J|£(r)| 2 -e(r)dr 3 

This expression is uniquely determined by the distribution of the dielectric constant e and the 
(unperturbed) electric field E and can, in most cases, be linearized in the small displacements 



u and reduced to a surface integral over the relevant resonator boundaries 1 10 20 1. 

In contrast to the unambiguous relation Q, the equivalent displacement x which is a-priori 
not a well-defined quantity, may be defined in an arbitrary way as discussed above. An impor- 
tant consequence of this ambiguity in the choice of x is the necessity to introduce an effective 



mass m e ff of the mechanical mode 1 10 21 1 which retains the correct relation for the energy U 



stored in the mechanical mode given a displacement x, 

U = l -m sn Q} m x 2 . (5) 

Evidently, rescaling x — > ax implies G — > G/a and m e ff — > m^ja 2 , so that the frequently 
quoted figures of merit G and m e ff are individually only of very limited use. It is possible to 
eliminate this ambiguity by referring to the vacuum optomechanical coupling rate 

go = G-x zvf (6) 



using terminology analogous to cavity quantum electrodynamics |22) . As expected from the 
Hamiltonian description Q, all optomechanical dynamics can be fully described in terms of 
this parameter. As a consequence, parameters such as the cooling rate or optical spring effect 
can be derived upon knowledge of go- For instance, the cooling rate in the resolved sideband 
regime is given by F = A{h ( ^jg 2 ) /K, while the power to reach the SQL scales with K 2 jg\ with 
the cavity loss rate K. Furthermore, the regime where go (n ) ^> fc,r m (T m is the mechanical 
damping rate) denotes the strong optomechanical coupling regime, in which parametric normal 
mode splitting occurs [3 23 24|. 

2. Determination of the vacuum optomechanical coupling rate 

Our intention here is to detail and verify an experimental technique to determine go in a simple 
and straightforward manner for an arbitrary optomechanical system. The technique we present 
relies on the fluctuation-dissipation theorem, from which the double-sided (symmetrized) spec- 
tral density of displacement fluctuations of a mechanical oscillator in contact with a thermal 
bath at temperature T can be derived [ 10 25|[26|, 



/ na \ 

S (O ..'Hi"— \ 2k B T J ... ' -' "'"<-->■ , 7) 

' 1 * - ' m eS - & + m eB (£22 - ^? m f + ' 



where ;« e ff is the effective mass corresponding to the choice of the normalization of x, as derived 
from ([5j, and the approximation is valid for large mechanical occupation numbers (n m ) m 
k#T /TiH m ^ 1. The displacement fluctuations are not directly experimentally accessible, but 
cause fluctuations of the cavity resonance frequency 0) c according to 



'For open systems, the integration boundaries should be chosen close to the surface of the resonator, to avoid 
divergence of the integral in the denominator. Alternatively, radiating fields escaping from the volume of interest have 
to be taken into account by a surface integral 1191. 



If such a spectrum can be measured, it is straightforward to integrate it over Fourier frequencies 
resulting in the interesting relation 

(5co c 2 ) = J^S aa (Q)^ = S (0(0 (a m ) I f =2(n m )g 2 , (9) 

which implies that the knowledge of the mechanical occupation number and the cavity reso- 
nance frequency fluctuation spectrum are sufficient to directly derive go. The former is usually 
a straightforward task (provided dynamical backaction [27 , 28 1 is ruled out, for example, by low 
probing powers), whereas the latter has to be implemented using a phase-sensitive measurement 
scheme as shown in figure [T] since the required changes of the optical resonance frequency 
change the phase of the internal and hence output field. 
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Fig. 1. Generic scheme for the determination of the vacuum optomechanical coupling rate. 
The cavity optomechanical system is probed with a monochromatic source of electromag- 
netic waves which is frequency-modulated by a known amount at a Fourier frequency close 
to the mechanical resonance frequency. A phase-sensitive detector generates a signal /, 
which is analyzed with a spectrum analyzer. As two possible examples, a homodyne de- 
tector measures the phase difference between the arms in an interferometer, one of which 
accommodates the optomechanical system (shutter S open). If the shutter S is closed, the 
transmission signal / of an optomechanical system is amplitude-dependent when the driv- 
ing wave is detuned with respect to the optical resonance. 

Various methods for such measurements are available, such as direct detection when the 
probing laser is detuned, frequency modulation |29|, polarization [30], or homodyne spec- 
troscopy pT[ . In general, cavity resonance frequency fluctuations are transduced into the signal 
in a frequency-dependent manner, such that the spectrum of the measured signal (e.g. a pho- 
tocurrent) / is given by 

Su{£i) = K(a)s vv (a) = / ^s a)0) (a), (io) 

where \j/ is the phase of the intracavity field, which gets modulated upon a change of the reso- 
nance frequency. Importantly, to exploit the relation ([9]), it is necessary to convert the spectrum 
Su(£l) back into a frequency noise spectrum, so the function K(Q.) = K(Q.,A, ri c , K,P) (where 
A = C0i — o c is the detuning from resonance, and tj c = k &x /k is the coupling parameter quan- 
tifying the ratio of the cavity loss rate due to input coupling Kg X compared to the total loss 
rate k) has to be calibrated in absolute terms. Experimentally, this can be accomplished by 
frequency-modulating the source that is used for measuring Sjj pT}[3"3~) . Under correct experi- 
mental conditions (cf. appendix) the phase modulation of the input field undergoes the exactly 
same transduction, so that 

s II (a)=K{a)Sn(a), (ii) 



where S^0(£2) is the spectrum of input field phase modulation. For the typical case of 
monochromatic input field phase modulation (E\ °< e~'^° cosamodt ), the equivalent spectral den- 
sity is given by S^(Cl) = 2n\ (5 (£2 - £2 mod ) + 8(Q. + Q. mod )) <j>%/2. 

A laboratory spectrum analyzer will, of course, not resolve this spectrum but convolute it 
with its bandwidth filter lineshape F(£l), for example a Gaussian centered at £2 = 0. If a time- 
domain trace is analyzed using a discrete Fourier transform, the filter function would be given 
by the squared modulus of the Fourier transform of the time-domain windowing function. We 
then have 

S™ eas (£2) = 2 • F(£i) * (K(a) (S vv {a)+S^(£l))) , (12) 

where a factor of 2 was introduced to account for the fact that most analyzers convention- 
ally display single-sided spectra (this factor is irrelevant for our final result, however). The 
filter function is usually normalized and can be characterized by the effective noise bandwidth 
(ENBW), for which 

d£2 



F(0)-ENBW= / F(£l)— = l, (13) 
J-™ 2% 

essentially absorbing differences in hardware or software implementations of spectral analysis 
(34). Note, however, that this description is still idealized. In reality, correction factors may 
apply depending on details of the signal processing (such as averaging of logarithmic noise 
power readings) in the respective analyzer used, see e.g. [35) . 

If the effective noise bandwidth is chosen much narrower than all spectral features in 
5^^(£2), the measured spectrum is approximately 

Sfr{&) « 2- A"(Q)S W (Q) + fK(a moi ) (FtQ-Omod) +F(n + £2 mod )) . (14) 

Assuming that the signal at £2 mo( i is dominated by laser modulation, the transduction function 
K can then be calibrated in absolute terms at this frequency, as in the relation 

sr s (±£w) « f m moi) F(o) = f . (is) 

all parameters except K(Q. mo d) are known. The spectral shape of K is then needed to derive the 
spectrum Syy of interest. This task can be accomplished either by repeating the measurement 
for all modulation frequencies of interest, or by theoretical analysis, relating it to the conditions 
under which the measurements are taken. 

As the simplest example, we may consider the case of direct detection, where the signal is 
the power P "transmitted" through the optomechanical system. For a given input power P ln , we 
find (see appendix) the transduction function: 

p lS ¥V PlS H ((£2 + A)2 + f c2/4)((£2-A) 2 + f c 2 /4)(A 2 + f c 2 /4) 2 ' k ' 

We assumed in this expression sufficiently low pump powers such that effects related to optome- 
chanically induced tranparency (OMIT) [36] are negligible. For the more advanced homodyne 
spectroscopy we calculate (see appendix), 

K H (0) ShH ShH 



4fc 2 77 2 £2 2 (£2 2 (l -2t] c ) 2 k- 2 /4 + (A 2 - (1 -2?7 c )k- 2 /4) 2 ) 
( (A - Q.) 2 + K 2 /A) ( (A + £2) 2 + K 2 /4) (A 2 + (l - 2r] c ) 2 K 2 /A) (A 2 + K 2 /4) ' 



(17) 



Advantageously, this method can also be applied with the laser tuned to resonance (A = 0), 
ruling out effects of dynamical backaction |27, 28 1 even for larger laser powers. In this special 
case 



Kh{£1) -WTW/4- (18) 

Another interesting case is A = ±i2 m in the resolved sideband regime [33 1 with i2 m 3> K, where 



MQ m )^4 JIc 2 (l + ^^]. (19) 



g 2 ( 16T7 c - 13; 

We would finally like to emphasize that in many experimental situations, the transduction 
function K(£l) is sufficiently constant over the range of frequencies, in which the mechanical 
response in (|9]l is significant (several mechanical linewidths around £l m ). In this case, it can be 
simply approximated to be flat, K(Q.) ps K(Q. m ) pa ^T(£2 mo d), and we get from equations (|9]), 
([10), ((TjE) and <[J3J> 

g2 ^ i 0o 2 QLd ^ eas (^ m )-r m /4 

^°~2(« m ) 2 ^(iWj-ENBW 



The last fraction in expression ( 20 1 can be interpreted as the ratio of the integrals of the dis- 
placement thermal noise and the calibration peak as broadened by the filter function. Note also 
that according to its definition the ENB W is given in direct frequency and not angular frequency 
as is r m . 

3. Application of the method 

As two examples, we have applied the discussed method to two different optomechanical sys- 
tems shown in figure [2] We first consider the case of a doubly-clamped strained silicon nitride 
beam (30 x 0.7 x 0.1 /im 3 ) placed in the near-field of a silica toroidal resonator (371. This 
scheme allows sensitive transduction of the nanomechanical motion with an imprecision below 
the level of the standard quantum limit 1 38 1. 

Note that the coupling rate in such a configuration depends on the strength of the evanescent 
field at the location of the nanomechanical oscillator, which increases exponentially if the os- 
cillator is approached to the toroid. The data in figure[2]show the single-sided cavity frequency 
noise spectral density S™y as (£l) = 5™ as (i2) • (H/lny- /K{Q) recorded using the homodyne 
technique. The spectrum was calibrated by laser frequency modulation at the Fourier frequency 
of 8 MHz. The spectrum needs not be processed further, since the transduction function for 
homodyne detection with A = is suffciently flat around the mechanical peak, r m <C K. How- 
ever, the cavity frequency noise spectrum as measured here contains also contributions from 
other noise mechanisms present in the silica toroidal oscillator, in particular thermorefractive 
noise [32| known to dominate frequency noise at lower (< 10 MHz) Fourier frequencies. The 
integration of the spectrum is therefore exclusively performed on the Lorentzian corresponding 
to the mechanical mode (gray area), which can be easily isolated by fitting the noise spec- 
trum with a model accommodating all known noise mechanisms including thermorefractive 



noise 1 38 39 . The result of this evaluation is (8(0%) = (2n ■ 530kHz) 2 , and the corresponding 
vacuum coupling rate is given by go = 2n ■ 420 Hz. 

Silica WGM microresonators such as spheres or toroids actually support structure-inherent 
micromechanical modes themselves [39, 40]. Figure [2] also shows an equivalent cavity fre- 
quency noise spectrum measured on such a cavity using direct detection. Again we neglect the 
frequency-dependence of K due to the very narrow mechanical peak (r m -C K,Q. m ). For the 
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Fig. 2. Measured (single-sided) frequency noise spectra S™ as (£l) (blue lines, v = ta c /27T) 
induced by a nanomechanical beam oscillator (8.3 MHz) coupled to a silica microtoroid 
|38| (left) and a radial breathing mode of a silica microtoroid (right), calibrated using laser 
frequency modulation at frequencies of 8 and 7 1.3 8 MHz, respectively. Orange lines are 
Lorentzian fits, the gray area underneath them yield the optomechanical coupling rates 
multiplied with the occupation number of the oscillator. 



evaluation of go, we remove backgrounds due to thermal noise from other modes, thermore- 
fractive fluctuations or shot noise, by fitting a Lorentzian model with a flat background term. 



Only the peak amplitude of the Lorentzian is then used with the relation ( 20 1. A vacuum cou- 
pling rate of go — 2n ■ 570 Hz is obtained, which is interestingly very similar to the vacuum 
coupling rate of the nanomechanical beam. 

4. Conclusion 

In summary we have presented a simple method which allows direct measurements of the vac- 
uum coupling rate go in optomechanical systems using frequency-modulation. We have ap- 
plied the method to two different devices, namely doubly-clamped strained silicon nitride beam 
(go/2n ss 10 3 Hz) and the radial breathing mode in a silica toroidal resonator (go/2n ~ 10 3 Hz). 
It is interesting to note how widely this coupling rate can vary in different optomechanical 



systems, for example, we estimate go/27!: ~ 10 2 Hz for (present) microwave systems |15|, 
go/2K ~ 10° Hz in lever-based optical systems fl2"||24) , and up to go 1 2% ~ 10 5 Hz in integrated 
optomechanical crystals [17]. Together with the optical and mechanical dissipation rates and 
resonance frequencies (which do also vary widely), the coupling rate provides a convenient 
way to assess the potential of optomechanical platforms for the different phenomena and ex- 
periments within the setting of cavity opto mechanics. 
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A. Transduction coefficients 



In this Appendix we show by direct calculations that the amplitude transduction coefficients of 
the input phase modulation and of the mechanical motion are the same both in the case of either 
phase homodyne or direct amplitude detection. 

The canonical equations describing dynamical behavior of the optomechanical system are 
the following: 

a = (iA-iGx(t) - K/2)a(t) + y/fycs^t), (21) 
x(f)+r m i(f)+£2^x(f) = -HG\a\ 2 , 

where a is the internal field in the cavity normalized so that \a\ 2 is the number of photons in 
the mode, i; n is the input field related to the input power as P in = fta) |sin| 2 , A is the detuning 
of the pump from optical resonance, K and T m are total optical and mechanical damping rates, 
?7 C is the coupling strength equal to the ratio of coupling and total optical losses and il m is 
mechanical resonance frequency. We neglect in the calculations below the effects of dynami- 
cal and quantum backaction of the optical probe field on the mechanical subsystem, implying 
sufficiently low input power. 

Assume that the mechanical degree of freedom of the cavity harmonically oscillates so that 

x(t) =x cos(i2 m f). (22) 

These oscillations will produce a periodic change of the optical eigenfrequency and hence 
the phase of the internal field with \j/Q = xqG / £l m . If these oscillations are small (|y/b| <C 1), the 
spectral components of the internal fields in the cavity and output field in the coupler will be: 

ax = w*fc*sf (o) ( i - ^•f+'H -^ - m^(-^) e+i g mt \ (23) 



JS?(Q) 



-i(A+Sl)+K/2' 



In the same way if the input field is weakly phase modulated (Sm °< e 'toa^^mcd')) to b tam a 
calibration signal: 

s $M = s iB (l - ^e-«W - ^ e +*w) , (24) 

i^S? (+flmod) „-iQ. m „ A t ^O-^(-^mod) +,o , 



a§ = W7cK"(^(0) 

s (j),out = s ij>,m — V ^ a ij) ■ (25) 

Dynamically these two modulations are not equivalent, but surprisingly produce the same trans- 
duction coefficient in terms of spectral densities. 

We discuss here two different methods to probe the output field, either by directly looking 
at optical power modulation on a detector with appropriate detuning A or by measuring phase 
oscillations using homodyne detection at zero detuning, with the local oscillator signal obtained 
by splitting the input power. From these relations it is possible to calculate the spectral densities 
of the power detected directly at the output of the detector, and the homodyne signal, once the 
spectral densities of mechanical (S XX (Q.)) or laser phase (S^(Q)) are known. 



A.l. Direct detection 

Directly calculating the modulation amplitude of the output optical power on a detector 



Pout(^) = ftcoo|soutln (omitting detector efficiency to convert to a photocurrent) we find that 



SppjO.) Spp(Q.) 
P?S vv (Cl) P*S H (Cl) 

477 c 2 K- 2 A 2 n 2 (i2 2 + K 2 (l - 77c) 2 ) 

~ ((i2 + A) 2 + K- 2 /4)((i2-A) 2 + K- 2 /4)(A 2 + )C 2 /4) 2 ' ( } 

In resolved sideband regime when Q. > \JlK the signal at frequency Q. is maximized when 
the pump is either tuned on the slope or on one of mechanical sidebands: 



A = ±^2tt 2 -K 2 -2 V / n4-2^~±|(\ + ^, 
A = ±iv // 2ii 2 -K- 2 + 2 v / i2 4 -2)C 2 i2 2 ~±i2^1-^^ , 



(27) 



producing the same signal 



r p (fl)=4T ?c :(l+ * (1 n2 l fc )) ). (28, 



(29) 



When Q. < V2k, the output is maximized when 

A = ±Wl2n 2 + 3K- 2 , 
6 

27fl 2 77 c 2 )c 2 (n 2 + K 2 (l~77 c ) 2 ) 

^° (£2) (ii 2 + f C 2 )3 • 

A.2. Balanced homodyne detection 

In this case, the signal on the detector is proportional to the field amplitudes slo and s; n falling 
on the output beam splitter, 

H . i 1 2 i 1 2 

^" l ^ out "^ Lo1 H 7f LO ~7f out| 

= i(' y LO s out ^ ^LO^out) = i( J LO*to - s LO s m) + iy/r^K (as£ - a*S LO ) 

= 2|jLo|kin| sin(0Lo - 0in) + J'x/tJc^^LO ~ «*^Lo) (30) 

It is important to note here that if the homodyne phase difference 0lo — 0in = then this 
method directly probes the phase of the internal field and hence the mechanical oscillations. 
This regime can be provided in case of A = 0. If the detuning is not equal to zero, the phase 
difference in an experiment is usually locked to have zero d.c. voltage. The local oscillator field 
is obtained by splitting the input power and is equal to: 

^,lo = — W ^°=*Loe^°, (31) 

(32) 



■s>,lo = s>,ine = Shoe \ t ~ e moa - —e 

Sir, \ Z Z 



where 0lo is chosen to suppress any d.c. voltage, 



In this case, again 

Shh Shh 



and 



(34) 



4>c 2 77 c 2 ^ 2 (n 2 (l - 277 c ) 2 >c 2 /4 + (A 2 - (1 - 2ifc) k 2 /4) 2 ) 
H[ ' ((A-i2) 2 + K- 2 /4)((A + i2) 2 + )C 2 /4)(A 2 + (l-277 c ) 2 K: 2 /4)(A 2 + K- 2 /4)' P ' 

If A = which appears the most favorable case for homodyne detection, then: 

i677 2 a 2 

^ (0) -^ + J c 2 /4- (36) 
Another interesting case is A = ±i2 m in case of resolved sideband regime. In this case 

M^)*4,?( 1+ ^-' 3 > ). ,37, 



